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The quenched hadron spectrum in the continuum obtained with the Wilson quark action in recent simulations 
on the CP-PACS is presented. Results for the light quark masses and the QCD scale parameter are reported. 
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1. Introduction 

Ever since the pioneering work of Hamber and 
Parisi[0 and of Weingarten||], the calculation of 
the light hadron spectrum has been regarded as a 
fundamental subject in lattice QCD. The exten- 
sive simulations on the GF11 computer^ have es- 
tablished that the quenched spectrum agrees with 
experiment within 5-10%. The next, and final, 
step within quenched QCD would be to further 
advance the calculational accuracy so that pos- 
sible deviations from the experimental spectrum 
are made manifest. In this article we report the 
main results recently obtained with the CP-PACS 
computer towards this goal. 

2. Simulation 

Our calculation is made with the plaquette 
gluon action and the Wilson quark action. We 
employ a spatial size of L s a rs 3 fm to avoid finite- 
size effects, and a -1 > 2 GeV to facilitate the con- 
tinuum extrapolation. These requirements lead 
to the selection of f3 and lattice sizes listed in 
Table [j]. Gauge configurations are generated 
with the over-relaxation and heat bath algorithms 
mixed in a 4:1 ratio. Hadron propagators are 
calculated for the five quark masses correspond- 



Table 1 

Parameters of simulation. 



p 


lattice 


a" 1 [GeV] 


L s a[fm] 


A con f 


5.90 


32 3 x 56 


1.97(1) 


3.21(2) 


800 


6.10 


40 3 x 70 


2.60(1) 


3.04(2) 


600 


6.25 


48 3 x 84 


3.12(2) 


3.03(2) 


420 


6.47 


64 3 x 112 


4.16(4) 


3.03(3) 


91 



*Talk presented by K. Kanaya at Lattice97, Edinburgh, 
Scotland, 22-26 July 1997. 



ing to m^/mp 0.75, 0.7, 0.6, 0.5 and 0.4, 
the last point being a step closer to the chiral 
limit than hitherto attempted. All unequal quark 
mass combinations allowed for degenerate u and 
d quarks are taken. We smear the quark source 
with the function \&(r) = A exp[— B(m p , m q ) r], 
with B w 0.33 fm. Hadron masses are extracted 
by uncorrelated fits together with a single elimi- 
nation jackknife procedure for estimating errors. 

3. Meson spectrum 

We find m 2 PS for pseudoscalar and my for vec- 
tor mesons to be well described by a linear func- 
tion of 1/K, which we adopt for our chiral fits. 
The physical point for the degenerate u and d 
quarks is fixed by 772^(135) and to p (769), and we 
use ?7iif (498) or 771^(1019) for the s quark. 

In Fig. we show the continuum extrapolation 
of strange vector meson masses with the vhk as 
input. A linear fit in a yields the continuum value 
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Figure 1. Continuum extrapolation of K* and <f> 
masses with the mjf-input. GF11 estimate for 
infinite volume S is plotted left ata= 0. 
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Figure 2. Meson hyperfine splitting. 

which is smaller than the experimental value by 
3% for K* and by 4% for cj), to be compared with 
the statistical error of 0.5%. If we take as the 
input for the s quark, to^* agrees with experi- 
ment to 0.6%, while mjc is higher by 9%. 

The origin of the discrepancy can be traced to 
the hyperfine splitting shown in Fig. ^. Our re- 
sults, which scale well, decrease too fast as meson 
mass increases. The slope b of the curve is related 
to the J parameter through J — (1 + b)/2. A fast 
decrease translates into a small value of J; we 
obtain J = 0.377(11) in the continuum limit. 

4. Baryon spectrum 

We show our typical chiral extrapolation for 
baryon masses in Fig. ||. We find all baryon 
masses to be linear in 1/K except for rnjy and 




Figure 3. Baryon masses as a function of 1/K 
for the light quark at (3 — 6.1. For non- 
degenerate baryons, results with the heavier s 
quark (m-n- /m p «0.75) are plotted. 
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Figure 4. Continuum extrapolation of baryon 
masses. Open symbols show GF11 results. 



toa, which show a clear negative curvature. We 
choose a fit cubic in 1/K for ttin, quadratic for 
m\, and linear for the other baryon masses. 

The negative curvature significantly lowers the 
nucleon mass compared to previous results at a 
finite lattice spacing, as shown in Fig. ^. The 
continuum value obtained from a linear extrapo- 
lation is smaller than experiment by 2.3%, albeit 
consistent within the statistical error of 3%. The 
mass of A is 6% higher, with a statistical signifi- 
cance of 1.7 standard deviations. 

The nucleon mass at each (3 depends on the 
form adopted for the chiral extrapolation. In- 
cluding the form from quenched chiral perturba- 
tion theory, the continuum values, however, agree 
within 3%. 



3 





1.6 




1.4 


> 






1.2 


CD 


E 


1.0 




0.8 




0.6 




0.4 



fl" 



I ■ • 



Experiment 

□ GF11 
• K-input 
<)>— input 



> 

<D 

> 4.5 

> 
CD 

a 



• Perturbative 






a wi local m 






. WI local exp(mj-1 






> Wl local 1-exp(-mJ 






■ Wl extended 








~ ~ - 


^ ^ ~ - - - 








0.0 0.1 0.2 


0.3 


0.4 



Figure 5. Quenched hadron spectrum in the con- 
tinuum limit. 



Figure 6. m= (ra u + md)/2 at 2 GeV in the MS 
scheme calculated by various definitions. 



The extraction of the continuum value is most 
uncertain for the nucleon, A, and X* masses, for 
which the errors in the continuum limit are about 
3%. For all the other channels a linear chiral ex- 
trapolation (quadratic for m\) followed by a lin- 
ear continuum extrapolation fits our data very 
well, resulting in errors of 1-2%. 

Our final result for the baryon spectrum is 
shown in Fig. ^ together with the meson spec- 
trum for completeness. 

In the octet the strange baryon masses are 
lower than experiment, by 5-8% using ttik as in- 
put, and by 3-5% for as input. The Gell- 
Mann-Okubo (GMO) relation, however, is well 
satisfied, at 1% for both cases. 

For the decuplet the GMO relation takes the 
form of an equal spacing rule. This is also well 
satisfied with our results, the three spacings mu- 
tually agreeing within 5-10%. The average spac- 
ings, however, are too small by 30% with ttik as 
input and by 20% for as input. 

The agreement of the baryon spectrum with 
experiment is generally better if we use as 
input, particularly for the decuplet. However, we 
have to remember that the K meson mass differs 
from experiment by 9% in this case, which is a 7 
standard deviation effect. 

5. Light quark masses 

We calculate light quark masses with the per- 
turbative definition m q a = Z m (l/K — l/K c )/2 
and with those based on the axial Ward iden- 



tity V pAy, = 2m q aP. The two methods give re- 
sults considerably different at finite lattice spac- 
ings. However, extrapolated linearly in the lat- 
tice spacing, they converge to a common value 

within 5% as shown in Fig. ^| For average u,d 

quark mass at 2 GeV in the MS scheme, we ob- 
tain m = 4.1(2) MeV. For the s quark mass, we 
find m s = 135(7)MeV and lll(4)MeV using to 
and rriK- The discrepancy reflects the quenching 
error in the strange meson spectrum. 

6. QCD scale parameter 

We estimate A^r using a P (3.40/a) and 
a Ms( 7r / a ) mcm ding the 3- loop correction in the 
A parameter and the scale determined from m p . 
The couplings are estimated from the average pla- 
quette with 2-loop perturbation theory. A linear 
continuum extrapolation of the two results con- 
verge and yields Aj^t = 230(5) MeV where the 
error includes extrapolation uncertainties. 
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